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ABSTRACT 

We study the origin of the extragalactic diffuse gamma-ray background using the data from the 
Fermi telescope. To estimate the background level, we count photons at high Galactic latitudes 
|6| > 60°. Subtracting photons associated to known sources and the residual cosmic ray and Galactic 
diffuse backgrounds, we estimate the Extragalactic Gamma-ray Background (EGB) flux. We find 
that the spectrum of EGB in the very-high-energy (VHE) band above 30 GcV follows the stacked 
spectrum of BL Lacs. LAT data reveal the positive (1 -|- z)^ , 1 < A: < 4 cosmological evolution of the 
BL Lac source population consistent with that of their parent population, FR I radio galaxies. We 
show that EGB at i? > 30 GcV could be completely explained by emission from unresolved BL Lacs 
if fc ~ 3. 

Subject headings: Gamma rays: diffuse background - BL Lacertae objects: general 



1. INTRODUCTION 

Study of diffuse background emission produced by 
faint sources with flux levels below the sensitivity of 
a telescope is commonly used to constrain the nature 
of source populations in the Universe and their cosmo- 
logical evolution. In the high-energy 7-ray band (HE. 
0.1-100 GeV) the diffuse extragalactic 7-ray background 
(EGB) was detecte d for the first time by SAS-2 satellite 
iFichtel et al.lll978[ ). was further s tudied by EGRET tele- 
scope on board oi CGRO mission (jSreekumar et al.lll998l : 
iStrong et al.ll2004 ) and, most recently, by the Large Area 
Telesc ope (LAT) on board of Fermi satellite (jAbdo et al.l 
l2010aD . 

It is often assumed that the dominant con- 
tribution to the EGB is given by distant Ac- 
tive Galactic Nuclei (AGN), in particular, 



tive Uaiactic iNuclei (AUIN), m p articu lar, blazars 
(iPadovan i et all Il993t iStecker et aTTll993t iCiang et al 
1995; S tecker fc Sa lamonI 119961 iMukherjee fc Chiang 



19991 1 iMiicke fc Poh]ll200d llnoue fc Totanj 120090 . How 



ever, a recent study by Fermi collaboration reveals that 
blazars might contribu te only a relatively small fraction 
of the HE EGB level (jAbdo et al.ll2010b[) . while a sig- 
nificant part of the EGB should be either explained by 
a yet unknown sourc e population or have a truly diffuse 
nature (see, however. IStecker fc Vente rs (2010,)). 

The EGB in the Very-High-Energy (VHE, 7-rays in 
the E > 100 GeV) range has never been measured. On 
one hand, the effective collection area of previous space- 
based 7-ray telescopes was not sufficient to achieve sig- 
nificant photon statistics in this energy band. On the 
other hand, the efficiency of cosmic ray background re- 
jection in the ground-based Cherenkov 7-ray telescopes, 
like HESS, MAGIC and VERITAS is not suflncient for 
detection of the isotropic diffuse EGB on top of the cos- 
mic ray background. Thus, the properties and the origin 
of VHE EGB remain largely unconstrained up to now. 



It is clear that the VHE EGB should contain a con- 
tribution from the unresolved point sources. The main 
candidate source class is, as in the case of HE EGB, 
that of blazars. At the same time, the VHE EGB could 
contain, apart from the contribution from unresolved 
extragalactic point sources, genuine diffuse components 
which could be produced via several mechanisms. For 
example, if the spectra of a large number of 7-ray-loud 
AGNs extend to the energies above 300 GeV, all the 
power emitted initially in 7-rays with energies higher 
than ~ 300 GeV is absorbed in the pair production of 
7-rays on the cosmological infrared a nd/or microwave 
backgrounds (jGould fc Schredeil Il967f ). Secondary in- 
verse Compton emission of electron-positron pairs de- 
posited in the intergalactic space in result of the pair pro- 
duction leads to generation o f diffuse extragalactic emis - 
sion in the VHE energy band ([Coppi fc Aharonianll997|) . 
Another mechanism which can lead to the generation of 
diffuse component of VHE EGB is electromagnetic cas- 
cade initiated in the intergalactic space by ultra-high 
energy cosmic rays (UHECR) interactin g with cosmic 
microwave background photons [Berezinskv fc Smirnovl 
[l97l[KaIa"shev et al.ll20"ol iBerezinskv et al.ll2011[ ). The 
cascade channels the power from the highest energies of 
about 10^" eV down to the ^0.1 GeV band in which the 
mean free path of the 7-rays becomes comparable to the 
size of the visible part of the Universe. 

Apart from the "guaranteed" (but, possibly, very 
weak) diffuse contributions, isotropic VHE 7-ray back- 
ground might contain contributions from "exotic" dif- 
fuse sources, like diffuse emission from annihilation of 
Dark Matter particles in the outer halo of the Milky Way 
galaxy and the annihilation signal accumulated from the 
dark matter halo s of all galaxies in the course of cosmo- 
logical evolution (jAbdo et al.ll2010cf ). 

Whatever are the sources of VHE EGB, they are scat- 
tered across the Universe, so that a significant con- 
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tribution to the flux is produced at redshifts z ~ 
1. A known effect of absorption of VHE 7-rays due 
to the interactions with infrared/optical Extragalac- 
tic Background Light (EBL) should lead to attenu- 
ation of the > 50 GeV signal produced by the 
sourc es at large red s hifts z ^ 1 (I Gould fc Schrederl 
1967t iKneiske et al I [20041: iFranceschini et alj 120081: 
Stecker et all 120061: iGilmore et all I2OO90 . This should 
leave an "imprint" on the VHE EBL spectrum, which 
should have the form of a gradual suppression with the 
increasing photon energy. Detecting a EBL suppression 
feature in the EGB spectrum would provide an impor- 
tant constraint on the (largely uncertain) evolution of the 
EBL density and spectrum up to redshifts z ~ 1. Such a 
constraint is otherwise difficult to obtain from the studies 
of individual extragalactic VHE 7-ray sources because of 
the limited signal statistics at the highest energies, espe- 
cially for the sources at significant redshifts. 

In what follows we discuss the measurement of the 
EGB derived from the data of Fermi/LAT telescope 
(| At wood et al.|[2009l ) . The measurement is obtained from 
the counting of photons at high Galactic latitudes, af- 
ter subtraction of the Galactic diffuse emission and the 
residual cosmic-ray background not rejected by the LAT 
data analysis software. We compare the measurement of 
EGB obtained in this way with the measurement previ- 
ou sly derived from the likelihood analysis of all-sky data 
bv lAbdo et all (l2010ah . 

The EGB flux above 30 GeV turns out to be com- 
parable to the flux in extragalactic VHE 7-ray sources 
resolved by Fermi. We find that the spectrum of EGB in 
this energy range follows the cumulative spectrum of the 
resolved sources. Dominant population of extragalactic 
VHE 7-ray sources is BL Lacs. Noticing the similarity of 
the spectrum of VHE EGB and of the cumulative BL Lac 
VHE 7-ray spectrum, we put forward a hypothesis that 
the VHE EBL is produced by unresolved BL Lacs with 
fluxes below the sensitivity of LAT. We explore this hy- 
pothesis and show that it could be valid if BL Lacs follow 
a positive cosmological evolution pattern, characteristic 
for other types of AGN, in particular for the parent pop- 
ulation of BL Lacs objects, Fanaroff-Riley type I (FR I) 
radio galaxies. 

2. DATA SELECTION AND DATA ANALYSIS 

For our analysis we consider all publicly available LAT 
data from August 4, 2008 to January 23, 2011. We pro- 
cess the data using Fermi Science ToolfQ. Wc filter the 
entire data set with gtselect and gtmktime tools following 
the recommendations of Fermi teanfl and retain only 
events belonging to "ultraclean" (P7ULTRACLEAN_V6) 
event class, which has minimal residual cosmic ray con- 
tamination. 

Estimate of the contribution of point sources to the to- 
tal flux requires separation of the photons coming from 
the point sources from those produced by the diffuse 
emission. Such separation is most straightforward for 
the photons with narrow point-spread- function (PSF). 
Taking this into account, select two sub-classes (which 
provide dominant contribution to the ultraclean events) 
with the most compact PSF, the sub-classes selected by 

^ |http : //fermi .gsfc .nasa.gov/ssc/data/analysis/sc ltools/ 1 
^ http: //fermi .gsfc. nasa.gov/ssc/data/analysls/scitools/ 



imposing the selection criterium EVENT_CLASS=65311 or 
32543. Other sub-classes of the ultraclean events have 
worse PSF. Point source contribution in these photons 
suffers from an additional uncertainty. Taking this into 
account, we restrict our attention to the subset of the 
ultraclean events with the best PSF. 

We retain events with Earth zenith angle 9^ < 100°. 
To estimate the flux from the photon counts we use gt- 
exposure tool. We consider only events at high Galactic 
latitudes, in the regions |6| > 60°. 

Our analysis is based on the so-called 
"Pass 7" selection of the L AT data (see 
http: //fermi .gsfc .nasa.gov/ssc/data/access/). 
However, use use a comparison of the Pass 7 data 
with the previous Pass 6 data selection in the estimate 
residual cosmic ray contamination of the set of events 
chosen for the analysis. The residual cosmic ray fraction 
in the P ass 6 data was studied in details by lAbdo et"al] 
()2010al) . Re-calculation of the residual cosmic ray 
fraction for any new selection of events, including the 
one considered in our analysis, could be done in a 
straightforward way as explained in Section 13.31 

3. DIFFUSE 7-RAY BACKGROUND 

Signal detected by LAT at high Galactic latitudes con- 
tains four types of contributions: emission from point 
sources, diffuse 7-ray emission from the Galaxy, EGB 
and residual cosmic ray background not rejected by the 
analysis software. To measure the EGB flux, one needs 
to separate the contributions from the four components 
in the overall signal in a given energy band. 

3.1. Point source contribution 

Point source component could be singled out in a 
straightforward way if the set of sources detectable in a 
given energy band is known. To define the set of sources 
we find the sources correlating with the arrival direc- 
tions of phot ons in each ener g y baii d, using the method 
described bv iNeronov et al. I (|2010[n . To calculate the 
total number of photons associated to the sources, we 
construct a cumulative distribution of photons as a func- 
tion of the distance 6 from the source and split it on the 
background and source contributions. The background 
contribution grows asymptotically as 0^ , while the source 
contribution asymptotically reaches constant. An ex- 
ample of the cumulative photon distribution around the 
source positions in the 12.5-25 GeV energy band is shown 
in Fig. m 

3.2. Galactic dijfuse emission contribution 

Contribution of the diffuse emission from the Galaxy 
should be found via a detailed fitting the all-sky photon 
distribution to an all-sky spatial and spectral template. 
This contribution is best constrained by the all-sky pho- 
ton distribution in the 0.1-10 GeV energy band, where 
event statistics is very high. Detailed fitting of the Galac- 
tic diffuse e mission to the data in the 0.1-10 GeV band 
was done bv lAbdo et al.l (j2010aD . In our analysis we rely 
on the best-fit model o f Galactic diffuse emission derived 
by lAbdo et al.l (|2010aD . This model is available in the 

^ See http: /www . Isdc .unige . ch/vhe/lndex . html for an up- 
dated version of the VHE source list. 
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Fig. 1. — Cumulative front photon distributions around point 
sources in 12.5-25 GeV energy band. Red points show the source 
photons, blue points show the background. Horizontal lines (from 
top to bottom) show the 100%, 95% and 68% levels. 

sky region of interest, |6| > 6 0°, see Fig. 6 of the Sup- 
plemental Material in the Ref. lAbdo et all ()2010aD . The 
uncertainties of t he Galactic d i ffuse em ission model are 
also discussed by lAbdo et"al] (j2010a[) . We take these 
uncertainties into account. 

The model consists of two main contributions: the 
"atomic hydrogen" component produced by interactions 
of cosmic rays with interstellar matter and "inverse 
Compton" (IC) component produced by inverse Comp- 
ton emission from cosmic ray electrons. Extrapolation of 
the atomic hydrogen component to the highest LAT ener- 
gies is straightforward: the pion decay spectrum follows 
the cosmic ray spectrum and extrapolation has the form 
of a simple powerlaw with photon index ~ 2.7, the same 
as the slope of the cosmic ray spectrum. This compo- 
nent gives a sub-dominant contribution above 100 GeV. 
Extrapolation of the IC component depends on the un- 
known shape of the (average over interstellar medium) 
cosmic ray electron spectrum at the energies above TeV . 
We have checked that in the model of lAbdo et al.l ()2010a|) 
the spectrum of IC component is consistent with the 
spectrum of IC scattering of the l ocal interstellar radia- 
tion field (iMoskalenko et aI]|2006D by electrons with the 
spectrum dNe/dE ~ E-^exp{-E/l TeV). This elec- 
tron spectrum consistent with t he cosmic ray electron 
spectrum observed on the Earth IjAckerman et al.ll2010l : 
I Aharonian et al.l 120081 ) . The IC spectrum produced by 
such electron population is shown by the cyan dotted 
line in Fig. The overall Galactic diffuse emission 
spectrum at high Galactic latitudes is then the sum of 
the atomic hydrogen and IC contributions, shown by the 
dashed black line in Fig. [2] 

The high-energy cut-off of the local cosmic ray elec- 
tron spectrum is most likely determined by the distance 
to the closest cosmic ray e lectron sources, e .g. to the clos- 
est pulsar wind nebulae (jAharonianl |2004[) . rather than 
by the intrinsic cut-off in the injection spectrum of elec- 
trons from the sources. This means that the local mea- 
surement of the high-energy cut-off of the cosmic electron 
spectrum does not provide a measurement of the high- 
energy cut-off in the injection spectrum of electrons. It 
is possible that Galactic cosmic electron sources inject 
electrons with energies much higher than ~ 1 TeV. This 
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Fig. 2. — Extrapolation of the spectrum of the Galactic diffuse 
emission at high Galactic latitudes |6t > 60° to the 100 GeV energy 
range. Blue and cyan dotted lines below 100 GeV show the con- 
tributions from cosmic ray interactions with interstellar medium 
an d inverse Compton scattering from cosmic electrons calculated 
bv lAbdo et al.l l|2010al ') . Continuation of the cyan dotted line above 
100 GeV is calculated assuming that inverse Compton emission is 
produced by electrons with a cut-off powerlaw spectrum with cut- 
off at 1 TeV. Grey dashed line is the sum of the cosmic ray and 
inverse Compton contributions. Solid grey line shows the overall 
diffuse emission spectrum in which the inverse Compton emission 
is produced by electron distribution without high-energy cut-off at 
1 TeV. Red thick solid line shows the spectrum used for subtrac- 
tion of Galactic component from the overall high Galactic latitude 
diffuse emission flux. 

possibility is shown by the solid grey line in Fig. [2]which 
shows the sum of the atomic hydrogen contribution with 
the IC emission from electrons without high-energy cut- 
off in the spectrum. The IC component still exhibits 
suppression at the energies 1 TeV because of the Klein- 
Nishina effect. 

To take into account the above mentioned uncer- 
tainty of the IC component we adopt an approximation 
dN-y/dE - E-^-^eyii>{-E/2 TeV) for the high Galac- 
tic latitude diffuse emission spectrum at > 100 GeV. 
This approximation is shown by the red thick solid line in 
Fig. [21 This spectrum lies exactly in the middle between 
the two extreme possibilities: TeV-scale high-energy cut- 
off in the cosmic electron spectrum and no cut-off in 
the cosmic electron spectrum. One should take into ac- 
count that the uncertainty of this approximation reaches 
~ 50% at the highest energies. We take this uncertainty 
into account in the calculation of the EGB spectrum, by 
adding it as a systematic error. The two extreme possi- 
bilities for the behavior of electron spectrum above 1 TeV 
(exponential cut-off exactly at 1 TeV and no cut-off at 
all) provide a good estimate of the overall uncertainty 
of electron spectrum in the interstellar medium in this 
range. The uncertainty of the shape of electron spec- 
trum dominates the uncertainty of the inverse Compton 
component of Galactic diffuse emission at high Galactic 
latitudes. 

3.3. Residual cosmic ray background contribution 

To estimate the residual cosmic ray background in 
the set of events selected for the analysis, we rely on 
the knowledge of residual The residual cosmic ray back- 
ground in the dataclean event class o f Pass 6 data 
is extensively discussed in lAbdo etlill ()2010aD . The 
residual cosmic ray background in the subset of Pass 7 
superclean events used in our analysis could be calcu- 
lated from the known residual cosmic ray background in 
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the Pass 6 dataclean events via a straightforard com- 
parison of statistics of events on- and off-point-sources 
in the two classes. 

First, the residual cosmic ray fraction in the Pass 6 
dataclean events should be calculated from the known 
suppression factor of cosmic ray event at transition from 
the diffus e event class to_the dataclean event class in 
Pass 6 (see lAbdo et"al] (|2010aD for the detailed discus- 
sion of the suppression factor). In each of the two event 
classes, the entire event set consists of a certain number 
of 7-ray events N^^i and a certain number of residual 
cosmic ray events, NcR,i where i stands for 3-1-4 or 4. 
Cleaning of the event set done to produce the dataclean 
event set from diffuse set results in rejection of a large 
fraction of the cosmic ray events, Ncr,a = ctCR^cR.a+i 
with acR <C 1. However, it results also in rejection of 
a number of true 7-ray events, so that = 
with a-y < 1. 

The suppression factor acR. is known function 
of energy from the Monte-Carlo simulations of cosmic 
ray and 7-ray indu ced events in the LAT detector by 
lAbdo et all (l20T0ah . The suppression factor could be 
found directly from the data set, by comparing statistics 
of events coming from the point sources in the diffuse 
and dataclean event classes (see section 13.11 above). 
In the calculation of all the photons associated to 
~ 10'^ poin t sources listed in the Fermi 2-year catalog 
(jAbdo et a l. 2010d) could be used. This provides very 
large event statistics so that uncertainty of is negli- 
gible. Knowing the total numbers of events in the two 
event classes Ntot,i one can resolve the system of equa- 
tions 

/ NcR,i/ acR + N^,4:/ = A^tot, 3+4 /-, n 

\ NcRA + N^A = NtotA 

with respect to Nj^^, Ncra to find the residual cosmic 
ray background in each energy bin for the dataclean 
event class. 

The residual cosmic ray fraction in the sub-class of the 
Pass 7 superclean events used in our analysis is then 
estimated in a similar way, once the residual cosmic ray 
fraction K4 in the point-source-subtracted set of the Pass 
6 events, Ncr.oSA = HiNotf^^ is known. 

Indeed, the transition from the Pass 6 dataclean 
events to the Pass 7 events belonging to the event classes 
65311 and 32543 leaves a fraction a^^^j of 7-ray events 
(actually, > 1 in a broad energy range around 10 

GeV) . It also suppresses or increases the residual cosmic 
ray background, so that the residual cosmic ray fraction 
in the off-source events changes from K4 to K7. The off- 
source events in the two classes are then the sum of the 
diffuse 7-ray emission photons and of the residual cosmic 
rays: 

f K;4A"off ,4 + A"7,off,4 = A"off ,4 

Knowing the statistics of the off-source events in the Pass 
6 and Pass 7 events, Nqq^^ and A^off,?, one could find the 
residual cosmic ray fraction in the Pass 7 data 

= 1 - Q-y.6->7(l - I^^I ^T^'^ (3) 
AoffJ 

The resulting estimates of the level of residual cosmic ray 



background for the events selected in the Pass 7 data in 
energy bins between 3 and 100 GcV are shown by the 
grey data points in Fig. [3] 

From Fig. |3]one could see that the contribution of the 
residual cosmic rays to the signal at 100 GeV is likely 
to be small. However, extrapolation of the estimate of 
efficiency of rejection of the residual cosmic ray back- 
ground much above 100 GeV is highly uncertain. It is 
possible that the efficiency of rejection of both the nu- 
clear and electron/positron component of the cosmic ray 
fiux drops because of the similarity of the cosmic ray and 
e^e~ pair tracks with large Lorentz factors. Inefficient 
rejection of the residual cosmic rays might lead to the 
contaminate the diffuse background signal and lead to a 
large over-estimation of the diffuse background fiux. Be- 
cause of this problem, we are able to only derive an upper 
limit on the EGB at the energies much above 100 GeV 
(for the energy band at 100 GeV we show a comparison 
between the 95% confidence level upper limit and the 
measurement). A proper measurement of the EGB flux 
at the highest energies accessible to LAT would require 
extensive Monte-Carlo simulations taking into account 
detector response (^ckermann 20l3)- 

3.4. Extragalactic ^-ray background spectrum 

EGB fiux could be found by subtracting the point 
source. Galactic diffuse and residual cosmic ray contri- 
butions to the overall number of events at high Galactic 
latitudes in each energy bin. The spectrum of EGB ob- 
tained in this way is shown by the red thick data points 
in Fig. 131 The error in the measurement of the EGB flux 
at the energies below 100 GeV has contributions from 
the uncertainty of the level of the residual cosmic ray 
background as well as from the systematic uncertainties 
of the Instrument Response Functions (IRF) and of the 
Galactic diffuse background in the relevant sky region. 
An additional contribution to the error in the VHE band 
is given by the statistical error arising from the low sig- 
nal statistics. Finally, one more uncertainty stems from 
the uncertainty of the shape of the cosmic ray electron 
spectrum in the TeV energy range, which propagates to 
the uncertainty of extrapolation of the inverse Compton 
component of the Galactic diffuse emission above 100 
GeV. 

In the same figure we compare the measurement of 
EGB obtained from the direct photon counting at high 
Galactic latitudes with th e results of the likeli hood anal- 
ysis of the all-sky data bv lAbdo et al.l ()2010al) . The two 
measurements agree well. 

Pink arrows at the energies above 100 Ge V show the 
uppe r limit on the VHE EGB derived by (jAckermannI 
[2010h using the likelihood analysis of the all-sky data. 
These upper limits agree with the upper limits derived 
from the direct photon counting (shown by the red arrows 
in Fig. 13 

As a matter of fact, the level of Galactic diffuse emis- 
sion at high Galactic latitudes turns out to be compa- 
rable to the level of EGB in the entire energy range 
> 10 GeV0 The Galactic diffuse emission contribu- 
tion to the total fiux at high Galactic latitude could not 

* There is no a-priori reason why the two fluxes should be nearly 
equal. Thus, the equality of the two contributions poses a "flne- 
tuning" problem which requires further investigation. 
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Fig. 3. — Estimate of the flux of isotropic component of 
diffuse emission obtained by the direct photon counting method 
(red thick line, data points and upper limits). For comparison, 
the spectrum of isotropic component of diffuse sk y emission, 
obtaine d using likelihood analysis at lower energies bv lAbdo et al.l 
Il2010al ). is shown as a red shaded r egion. Pink upper limits 
above 100 GeV are from [XckermannI I I2010I ). Black data points 
show the total point-source-subtracted flux from the North 
and South Galactic pole regions at \b\ > 60°. Solid line er- 
rorbars show statistical error. Dashed line errorbar show the 
systematic error at the level of ~ 20% stemming from the uncer- 
tainty of the the Instrument Response Functions (IRF) of LAT (see 
ihttp : //f erml . gsf c .nasa. gov/ssc/data/aiialysis/LAT_caveats .html ). 
clack horizontally shaded region shows the point source subtracted 
flux in the Galactic Pole regions found by lAbdo ct al.i i2010al ). 
Grey data points and grey curve show the estimate of the residual 
cosmic ray background in the event set used in this analysis. The 
res idual cosmic ray b ackground level in the data set considered 
bv lAbdo et al.l Il2010al ) is shown by the grey shaded region. Blue 
shading shows the Galactic diffuse emission in the North/Sout h 
Galactic Pole regions \b\ > 60° derived by lAbdo et all I l2010al '). 
Blue line shows the Galactic diffuse background spectrum. 

be negligibly small in the 100 GeV band. Taking this 
into account, it is not surprising that our estimate of the 
VHE EGB flux is somewhat lower than the total diffuse 
emission flux at high Galactic latitudes an d is, respec- 
tively, lower than the upper limit derived by lAckermannl 
(IMO) . 

It is useful to note that extrapolating the EGB spec- 
trum as a powerlaw spectrum to E > 100 GeV band 
would give the spectrum consistent with the data above 
100 GcV. With the current LAT exposure, there is still 
no evidence for suppression of VHE EGB flux due to ab- 
sorption on EBL. A larger exposure time is needed to 
verify the presence of the feature. As it is mentioned in 
the Introduction, suppression of the flux above 100 GeV 
due to the absorption of VHE 7-rays on the EBL is ex- 
pected if the EGB is accumulated over the cosmological 
distance scale. Detection of such suppression would be 
an important test of the origin of EGB. 

4. E>30 GEV EXTRAGALACTIC 7-ray BACKGROUND 
FROM POINT SOURCES 

As it is mentioned in the Introduction, different types 
of point and diffuse sources could contribute to the EGB 



Fig. 4. — Cumulative 7-ray flux from different classes of resolved 
point sources. Source class is marked to the left from each curve. 
Green shaded are a shows the point s ource flux at high Galactic 
latitudes found bv lAbdo etaHl Il2010al) . Red curve shows the EGB 
spectrum from Fig. [S] which includes only unresolved sources. 

in the VHE band. The main class of extragalactic point 
sources detected by Fermi is blazars, which are divided 
onto two sub-classed: BL Lac type objects and Flat 
Spectrum Radio Quasars (FSRQ). Over the first year 
of operation LAT has det ected some ^ 700 s uch objects 
above 100 MeV energy ()Abdo et al.l 120 lObl) . BL Lacs 
and FSRQ have somewhat different spectral character- 
istics in the 7-ray band, with the spectra of BL Lacs 
being syst ematically harder than the spectra of FSRQ 
(jAbdo et a l. 2010b). Hardness of the spectra of BL Lacs 
implies that they might produce significant contribution 
to the overall 7-ray flux in the VHE band. In fact, most 
of the extragalactic VHE 7-ray sources detected up to 
now by the ground based 7-ray telescopes sensitive above 
100 GeV are BL Lacfl 

Fig. [3] shows the breakdown of the point source contri- 
butions to the high Galactic latitude flux by the source 
type. One could clearly see that the dominant contribu- 
tion is given by BL Lac objects which provide > 90% of 
the total point source flux above 30 GeV. The cumula- 
tive spectrum of the other major blazar class, FSRQ has 
a high-energy cut-off at ^ 10 GcV so that FSRQ contri- 
bution to the point source flux is negligible in the VHE 
band. From Fig. [4] one could see that the total point 
source flux calculated from the cumulative photon distri- 
bution around stacked point sources in the high Galac- 
tic latitude regions (sec Section lO]) is in a good agree- 
ment with the total p oint source flux calcu lated using the 
likelihood analysis by lAbdo et al.l ()2010aD , shown by the 
green shaded region. Above 50 GeV 90% of source pho- 
tons come from BL Lacs and 10% from "Other" Fermi 
sources, which are dominated by not-identified sources 

^ For the catalogs of ext ra galactic VHE 7-ray 
sources see e.g. http : / /tevcat .uchicago . edu| and 

http : / /www . i sdc. unige . ch/ vhe / index . html 
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Fig. 5. — Number of detected VHE photons as a function of 
redshift in the 6.25-12.5 GeV (blue dotted histogram) 25-50 GeV 
(green dashed histogram) and 100-200 GeV (red solid histogram) 
bands. Also shown are the distributions of photons with the red- 
shift expected for different laws of BL Lac cosmological evolution. 

with some contribution from nearby AGN's. 

The EGB spectral shape above 30 GeV follows the cu- 
mulative point source spectrum. This observation leads 
to a conjecture that the VHE EGB is produced by al- 
ready known type of VHE 7-ray point sources with fluxes 
below the sensitivity of LAT. Since the dominant source 
class in the VHE band is that of BL Lacs, a more pre- 
cise conjecture is that the VHE EGB is produced by the 
unresolved BL Lacs. 

5. BL LAC CONTRIBUTION TO VHE EGB 

In the unification schemes of AGN BL Lac objects 
are identified with the FanarofF-Riley type I (FR I) ra- 
dio galaxies with jets aligned with the line of sight 
(jUrrv fc Padovanilll995l ). This implies that the cosmo- 
logical evolution of the I3L Lacs should follow that of the 
FR I radio galaxies. Recent studies of the cosmologi- 
cal evolution of FR I galaxies show that they experience 
"positive" cosmological evolution, which is usually de- 
scribed in terms of luminosity or comoving source density 
evolution as the increase of either average source lumi- 
nosity or the average comoving source density with the 
redshift (1 -I- z)'', k > 0. Different recent studies 

find somewhat different values of fc, depending on the 
analyzed radio galaxy samples and different assumptions 
about the evolution typ e (luminosity or density), with 
fc ran ging in 1 < fc < 3 (iSadler et a lJl2007t iHo dge et alJ 



l2009tlSmolcic et al.ll2009D . Since BL Lacs are just the FR 
I galaxies specially oriented with respect to the line of 
sight, their cosmological evolution follows the evolution 
of FR I galaxies, with the increasing source luminosity or 
spatial density with the redshift. This implies that sig- 
nificant flux should be produced by the sources at large 
redshifts, z ^ 1. It is possible that most of individual 
sources at large redshifts are too weak to be significantly 
detected by LAT, but collective emission from all the set 
of BL Lacs at high redshifts gives a significant contribu- 
tion to the EGB. 

Dependence of the total flux of the BL Lac popula- 
tion on the redshift could be found from the follow- 
ing straightforward calculation. Let us consider the 
total flux produced by sources at redshift z in a red- 
shift interval Az. This redshift interval corresponds to 



^/fl^+T^^,JY+z)^ is the expansion rate of the 
Universe fllled with matter and cosmological constant 
with today's densities flm and f2A- 

As an example, we take the case of "pure luminos- 
ity" evolution with the average source luminosity in- 
creasing as {1 + z)'' and conserved comoving source den- 
sity n{z) = uq = const. The number of sources in a 
spherical layer of thickness Az is ANg = Aim^r'^Ar. 
Each source produces the flux in a given energy band 
F ~ (1 -I- z)''~^ / (47rr^), where F is the photon index, the 
factor (1 -f- z)^~^ describes the change in the number of 
photons in a given energy band due to the cosmological 
redshift of the photon energies. One power of (1 + z) 
is compensated by the time delay between subsequent 
photons. 

The flux from the sources at large redshifts is affected 
by absorption of VHE photons on EBL. For example, at 
z ~ 1.5 the absorption modifies source spectrum above 
the energy i? ~ 50 GeV, if one assum es the EBL evo- 
lution calculated by iFranceschini et al] (j2008i) . Absorp- 
tion on EBL leads to suppression of the flux by a factor 
exp (— r(i?, z)) where t{E,z) is the optical depth with 
respect to the pair production. 

The overall flux from the sources in the redshift interval 
Az is 



AF{E,z) 
A'z 



FAN, 



(1 + ^) 



k-r„-T{E,z) 



y/nA + n,„ii + zf 



(4) 



Fig. [S] shows the number of 7-rays as a 
function of source redshift. For this we used 
BL Lacs with known red shifts from Veron&Veron 
(jVeron- Getty fc VeronI |201CI[ ) catalog complemented by 
BL Lacs detected by LAT, but not listed in the 
VeronfcVeron catalog. Only sources with |&| > 10° 
were considered. Here we plot photon distributions from 
Fermi BL Lacs in the three energy bands: 6.25 — 12.5 
GeV, 25 - 50 GeV and 100 - 200 GeV. One can see that 
at lower energies -E < 50 GeV a significant flux is pro- 
duced by BL Lacs at large redshifts up to z = 1.5. At the 
highest energies only contribution from nearby sources at 
z < 0.7 is present. Two effects might explain the deficit 
of high-redshift sources at high energies. First, the flux 
at the highest energies is suppressed by absorption on 
EBL. Next, the photon statistics in the highest energy 
bin is low so that sources contributing to the flux in the 
6.25-12.5 GeV bin produce less than one photon in the 
100-200 GeV bin. 

In the same figure we also show the expected de- 
pendence of the number of photons on the redshift ex- 
pected in different evolution models, Eq. The mod- 
els for cases fc = 1,2,3 are shown with magenta lines for 
6.25 — 12.5 GeV energy band. We normalize the mod- 
els to the number of photons in the first redshift bin, in 
which we have the most complete knowledge of the BL 
Lac population. 

From the comparison of the evolution models with the 
data one might get an impression that (1 + z) model is 
more consistent with the data than the models assum- 
ing faster evolution. However, the histogram on Fig. [5] 
does not take into account photons from BL Lacs with 
unknown redshifts. These BL Lacs produce about 30 % 
of all cumulative BL Lac flux. This means that at least 
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Fig. 6. — VHE EGB produced by unresolved BL Lacs under 
different assumptions about the cosmological evolution of BL Lac 
population (the evolution law is marlicd to the left of each curve). 
Red data points show the EGB spectrum from Fig. |3] 

30% contribution to the overall flux (integrated over all 
redshifts) is missing in Fig. [5] The model with evolu- 
tion k = 1 predicts the total number of photons which is 
~ 3(7 below the total number of photons in BL Lacs with 
known and unknown redshift together in the 6.25 — 12.5 
GeV energy band. In the energy band 3.125-6.25 GeV 
the undcr-prediction of the total number of photons from 
BL Lacs in the fc = 1 model is at > 5cr level, which means 
that the model is efficiently ruled out. Thus, Fermi LAT 
observations of BL Lac objects indicate that BL Lacs 
have positive cosmological evolution with > 1. 

In all other cases fc > 1, the discrepancy between the 
observed and expected number of photons from BL Lacs 
starts already at small redshifts, z > 0.2. The "missing 
BL Lac" 7-rays could come either form BL Lacs with 
unknown redshifts or from Fermi sources which are not 
yet identified as BL Lacs or, finally, from BL Lacs with 
fluxes below the sensitivity of LAT. 

Although individual high redshift BL Lacs would not 
be detectable by LAT, cumulative flux of these BL Lacs 
could give significant contribution to the EGB. Fig. IH] 
shows the contributions from "missing BL Lacs" at high 
redshifts up to 2; = 1 expected in four different models of 
cosmological evolution of BL Lac / FR I population. To 
calculate this contribution, we have normalized AF/Az 
distribution on the measured fiux of BL Lacs in the red- 
shift bin <z < 0.1. For the two last bins,100-200 GeV 
and 200-400 GeV, the statistics of the BL Lac signal is 
too low to normalize AF/Az by the fiux in the first red- 
shift bin. To estimate the normalization of AF/Az in 
these energy bins we have assumed that the average BL 
Lac spectrum extends as a powerlaw with photon index 
—2 up to the 400 GeV energy range. The resulting statis- 
tics of the signal in the last 200-400 bin in Fig. IHlis low 
and is subject to large fluctuations. 

From Fig. |6] one could see that the conjecture that 



EGB is produced by collective emission from distant BL 
Lacs is valid if BL Lacs experience positive luminosity 
or density evolution of the form (1 -|- z)*"' with fc ~ 3. 
Slower or faster cosmological evolution with fc = 2 or fc* = 
4, respectively, would under- or over-produce the EGB. 
Since evolution fc = 4 predicts number of photons more 
then diffuse gamma-ray background (see Fig. [5]) it is also 
excluded. Thus, the LAT data impose a constraint on the 
cosmological evolution of the BL Lac/ FR I population 

1 < fc < 4 (5) 

From Fig. |6] one could also see that if fc = 3, all the 
flux of EGB above 30 GeV could be explained by a cu- 
mulative emission from the BL Lac population. At the 
same time, if fc is significantly smaller than 3, as it is 
suggested by some recent studies o f the evolution of th e 
parent population of FR I galaxies (jSmolcic et al.|[2009l) , 
significant part of the VHE EGB fiux should come from 
a yet unknown source population or have a truly diffuse 
nature. 

7-ray flux from high-redshift BL Lacs is modified in the 
VHE band by the effect of absorption on EBL. The model 
spectra shown in Fig. IHtake this effe c t into account. We 
use the model of iFranceschini et al.l (|2008[ ) to estimate 
the attenuation of the VHE 7-ray flux from redshifts up 
to z ~ 1.5. Model curves shown in Fig. |6] assume that 
the average intrinsic spectrum of BL Lacs does not have 
a high-energy cut-off up to 400 GeV. The observed 
suppression of the flux above 100 GeV is explained only 
by the effect of absorption on EBL. 

6. DISCUSSION AND CONCLUSIONS 

In this paper we have derived a measurement of EGB 
in the 10-400 GeV energy range from the analysis of 
Fermi/LAT data in the North and South Galactic Poles 
regions, \b\ > 60°. Our approach was to count all 
the protons detected by LAT in this region and esti- 
mate the number of counts from the point sources, from 
the Galactic diffuse emission, form the residual cosmic 
ray background and from the EGB. Subtracting the 
source. Galactic diffuse and residual cosmic ray back- 
ground counts from the total number of counts in the 
North and South Galactic Pole regions we derived the 
EGB spectrum shown in Fig. [31 

Comparing the spectrum of EGB in the > 30 GeV 
energy band with the spectrum of cxtragalactic point 
sources in the same energy band (Fig. |4]) , we have noticed 
that the two spectra closely follow each other. Based on 
this observation, we have put forward a conjecture that 
the EGB above 30 GeV is explained by the unresolved 
BL Lacs, which give dominant contribution to the extra- 
galactic point source flux in this energy band. We have 
demonstrated that this conjecture is consistent with the 
EGB measurement provided that BL Lacs follow positive 
cosmological evolution with the overall power of emis- 
sion from the source population increasing as (1 -I- z)^ 
up to z 1 (Fig. [6]). Such cosmological evolution is 
roughly consistent with the measurements of cosmolog- 
ical evolution of FR I radio galaxies which are believed 
to be the parent population of BL Lac type objects and 
are also observed to have positive cosmological evolu- 
tion of the form (1 -I- z)'' with an uncertain value of fc 
between 1 and 3 jRigbv et al.|[200llSmolcic et~alll2009l: 
iSadler et al.ll2007D . At the same time, it is opposite to 
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the negative cos mological evolution of the high-energy- 
peaked BL Lacs ()Giomnii et aL|[r999ll200il ). which con- 
stitue a sub-class of the GeV-TeV 7-ray emitting BL Lacs 
considered in our analysis. 

If the VHE EGB is indeed produced by distant BL 
Lacs at rcdshifts up to z ^ 1, LAT will not be able to 
resolve it into point sources. Indeed, the brightest BL 
Lac on the sky, Mrk 421 produced ~ 30 photons above 
100 GeV. If the positive cosniological evolution of BL 
Lac population is mostly due to the increase of the co- 
moving source density rather than increase of the typical 
source luminosity, the brightest BL Lacs at redshift 2^1 
produce ~ 10~^ 7-rays in LAT over some 2.5 years of ex- 
posure. This means that LAT would not collect sufficient 
photon statistics to detect distant BL Lacs individually. 
If the (1 -I- zY evolution is mostly due to the increase of 
the average source luminosity, BL Lacs at redshift 1 have 
an order-of-magnitude higher luminosity than local BL 
Lacs. However, even with higher luminosity, they pro- 
duce on average 0.1 photon in LAT, so that they are still 
not individually detectable. 

If the real value of k is much below k = 3, as indi- 
cated by a recent study by iSmolcic et al.l ()2009l ) , emis- 
sion from the unresolved BL Lacs would not explain the 
VHE EGB flux and there should be another source class 
or a mechanism of production of diffuse emission which 
would account for the EGB. Such a mechanism might, 
in fact, be indirectly related to the BL Lac population. 
Most of the power output from BL Lacs at the ener- 
gies above 100 GeV is converted into the electromagnetic 
emission from 7-ray induced cascade in intergalactic 
medium. The intrinsic spectra of BL Lacs (and of FR I 
radio galaxies, such as M 87 and Gen A (jAharonian et al."] 
[200I lAharonian et al.l 12001 )) are known extend up to 
~ 10 TeV. Typical energy of the cascade photons 
which are produced via inverse Compton scattering of 
GMB photons by the e+e" pairs deposited in the in- 
tergalactic medium is ~ im[E^JlOTeVf GeV 
()Neronovfc Semikozl 120091 ). If the intrinsic source lumi- 
nosities in the 1-10 TeV range are comparable to the 
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luminosities in the 10-100 GeV, total flux of the cas- 
cade emission in the 10-100 GeV band is expected to be 
comparable to the point source flux, so that the cascade 
emission could give signific ant contribution to the EGB 
(jGoppi fc Aharonianlll997[ ). This is consistent with the 
observation that the VHE EGB level is comparable to 
the cumulative extragalactic point source flux in the 10- 
100 GeV band observed by LAT. 

The only possibility to test the hypothesis of BL Lac 
origin of EGB would be to use deep observations with 
ground-based 7-ray telescopes. Ground based 7-ray tele- 
scopes, which are sensitive in the VHE energy band, have 
much larger collection area above several hundreds of 
GeV and, as a consequence, could detect much weaker 
sources, than LAT. The flux from Mrk 421-like BL Lacs 
at the redshift z ~ 1 is ~ 10" ^'^ erg/cm^s. If the 
cosmological evolution of BL Lacs is due to increase 
of the average source luminosity with redshift, bright- 
est BL Lacs at redshift z ^ 1 might produce fluxes up 
to 10~^^ erg/cm^s at 100 GeV. At the energies around 
E < 100 GeV this flux is not strongly attenuated by the 
absorption on EBL. The energy E ~ 100 GeV is around 
or below the low energy threshold of the current genera- 
tion Cherenkov telescopes, like HESS, MAGIC and VER- 
ITAS. However, next g eneration facilities, Ch erenkov 
Telescope Array (CTA) (ICTA Consortium! [2OIOI) or 5@5 
(jAharonian et al.l 1200 ID are expected to have an energy 
threshold significantly below 100 GeV. Their sensitivity 
could be sufficient to resolve the VHE EGB into point 
sources, at least in the case when the cosmological evo- 
lution is mostly luminosity, rather than source density 
evolution. 
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